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ABSTRACT 

Absorption and emission lines originating from the C-H stretching manifold of gas phase CH4 were 
discovered in the high resolution (_R=25,000) infrared L band spectrum along the line of sight toward 
NGC 7538 : IRS9. These observations provide a diagnostic of the complex dynamics and chemistry 
in a massive star forming region. The line shapes resemble P Cygni profiles with the absorption and 
emission components shifted by ~ 7 km s""'^ with respect to the systemic velocity. Similar velocity 
components were observed in CO at 4.7 /im, but in contrast to CH4, the CO shows deep absorption 
due to a high velocity outflow as well as absorption at the systemic velocity due to the cold outer 
envelope. It is concluded that the gas phase CH4 abundance varies by an order of magnitude in this 
line of sight: it is low in the envelope and the outflow (X[CH4]< 0.4 x 10^®), and at least a factor of 
10 larger in the central core. The discovery of solid CH4 in independent ground and space based data 
sets shows that methane is nearly entirely frozen onto grains in the envelope. It thus appears that 
CH4 is formed by grain surface reactions, evaporates into the gas phase in the warm inner regions of 
protostellar cores and is efhciently destroyed in shocks related to outflows. 

Subject headings: infrared: ISM — ISM: molecules — ISM: abundances — stars: formation — stars: 
individual (NGC 7538: IRS9) — astrochemistry 



1. INTRODUCTION 

Among the suite of molecules observed in interstel- 
lar and circum-protostellar media*, methane is relatively 
poorly studied. CH4 has no permanent dipole mo- 
ment, and therefore cannot be observed by pure rota- 
tional transitions at radio wavelengths. Its strongest 
fundamental [v = 1 — 0) ro-vibrational transitions, i.e. 
the C-H stretching and the C-H bending modes, 
occur at 3.32 and 7.67 /J,m, respectively. Strong tel- 
luric absorption at these wavelengths hinders interstel- 
lar CH4 studies. Gas and solid phase CH4 have been 
positively identified at 7.67 /Ltm in three sight-lines: the 
massiv e protostars NGC 7388 : IRS9, W 33A, and GL 
7009 S llLa,cv et a,1. 1 99l!: IBoogert et a ,l. 1 996.T997'. [T99^ 
iDartois et alJll998 'l. CH4 was detected only in the solid 
state in low resolution Infrared Space Observatory (ISO) 
observations at 7.67 /im toward ~ 30 low mass protostars 
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OPERATED AS A SCIENTIFIC PARTNERSHIP AMONG THE 
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^ California Institute of Technology, Division of Chemistry and 
Chemical Engineering, Pasadena, CA 91125 

* see the compilation at 
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iWhite et ani2000t ICurtler et al.ll200l [Alexander et alJ 
I2003|) . theGalactic Center (Chiar et al. 2000^. and an 
external galactic nucleus ((Spoon et al. 2001) . A recent 
low resolution Spitzer Space Telescope spectrum shows 
absorption by solid CH4 toward the young, low mass pro- 
tostar IRAS 08242 -H5100 (Noricga-CrcsDO et al. 200| 
iBoQgert et a"l]l200^ . 

Observations of CH4 provide insight into the basic 
principles of astrochemical networks. In the presence 
of atomic C and H, CH4 is rapidly formed on cold 
grains , just as H2O is formed from at omic O and H 
(e.g. iBrown. C harnlev. fc Millar I Il988f) . Perhaps sur- 
prisingly, the observed gas and solid state CH4 abun- 
dances are low, not m ore than a few percent with re- 
spect to CO and H2O l|Lacv et alJ 119911: iBooeert et alJ 
Il998j) . This points to relatively low atomic C abun- 
dances at the time of CH4 formation, when most C 
is already locked up in CO. In this same scenario the 
high CH3OH abundances in several lines of sight may 
be explai ned by h ydrogcnation of abundantly p resent 
CO tP artois et al.l ll999: Pontoppida n et al J 12003) . and 
large CO2 abundances by oxidation of CO. The dom- 
inance of a CO- and 0-based, rather than C-based, 
chemistry is stressed by the profiles of solid CO2 ab- 
sorption bands which show an i ntimate mixture of CO2 , 
CH.-hOH, and H2O i n the ices ijGerakines et al. I Il999t 
iBoQgert et al. ll200nD . Nevertheless, even at relatively 
low abundances methane is proposed to be the starting 
point of a rich chemistry leading to the complex organic 
molecules observed in dense clou ds like the TMC-1 ridge 
iMarkwic k" Millar, fc CharnlevI [20001. To further test 
the role of CH4 in interstellar chemistry, basic diagnos- 
tics, such as gas phase abundances, gas/solid phase abun- 
dance ratios, gas phase temperatures and abundance gra- 
dients need to be measured in more lines of sight. 

In this work we report the discovery of absorption and 
emission features arising from the 3.32 /im C-H stretch- 



2 



Boogert, Blake, & Oberg 



ing mode of gas and solid phase CH4. A previous attempt 
to detect CH4 at this wavelength was unsuccessful and in- 
dicated an abund ance of < 10~^ with respect to H in the 
Orion BN object ijKnacke et all 198511 . We study the line 
of sight toward the massive protostar NGC 7538 : IRS9, 
because its high radial ve locity of Visr = —57 km s~^ (e.g. 
Ivan der Tak et al. 12000(1 facilitates the separation of tel- 
luric and interstellar lines. NGC 7538 : IRS9 is a well 
studied massive protostar, known for its large ice and gas 
columns l|Whittet et al. 1199 6: Mitchell ct al. 19901 and 
a prominent molecular outflow (Mitchell, Maillard, & 
Hasegawa 1991). The R(0) and R(2) lines of the bending 
mode of gas phas e CH4 and solid sta te CH4 were detected 
in this source bv iLacv et al.l l|199lD . Unambiguous con- 
firmations of these detections were obtained with space 
based observations (Booge rt et al...l99 6. 1998). Here we 
independently confirm the presence of CH4 by detecting 
many lines from the 3.32 fj,m stretching mode as well as 
the ice band. The high spectral resolution of these obser- 
vations permits a line profile analysis in which we deter- 
mine gas phase conditions in the various environments 
along the line sight, such as the envelope, the hot core, 
and the outflow. This, in turn, elucidates the formation 
and destruction pathways of CH4 in these different envi- 
ronments. 

The Keck/NIRSPEC and ISO/SWS observations are 
discussed in § |21 and presented in § 13.11 A proflle analy- 
sis of the gas phase absorption and emission CH4 lines, 
in conjunction with previously reported high resolution 
CO observations, is presented in SS I3.2.1l and 13.2.21 The 
ice band profile is discussed and compared with the C-H 
bending mode in § 13.31 The implications of these re- 
sults on the formation and evolution history of the CH4 
molecule in protostellar environments are discussed in 
§ 14.11 The origin of the low velocity, warm gas, respon- 
sible for the P Cygni-like line profiles, is further con- 
strained in § 14.21 We conclude with future prospects in 

m 



2. OBSERVATIONS 



2.1. Keck/NIRSPEC R^25, 000 Spectroscopy 

The high resolution L band spectrum of the massive 
protostar NGC 7538 : IRS9 was obtained w ith the NIR- 
SPEC spectrometer l|McLean et al. Ill998|) at the Keck 
II telescope atop Mauna Kea. Observations were made 
on five different nights with clear and dry weather con- 
ditions: UT 2002 December 17 and 18, 2003 July 10 and 
11, and 2003 November 2. NIRSPEC was used in the 
echelle mode with the 0.43x24" slit, providing a resolv- 
ing power of i? = A/AA = 25,000 (^12 km s'^). The 
2.9-4.1 fim region was almost fully covered in five echelle 
settings, each covering five Nyquist sampled wavelength 
orders of ~ 0.05 /xm on the 1024^ pixel array. Two small 
regions centered on 2.97 and 3.12 /xm were not covered. 
The on-source integration time per setting was 20 min- 
utes. Several CII4 ro- vibrational lines as well as the solid 
CII4 band at 3.32 /im were observed multiple times to 
verify their reahty. At Visr = -57 km s"^ for NGC 7538 
: IRS9, the Doppler velocity with respect to earth was 



—53.5, —82.4, and —64.0 km s^^ for the three observing 
runs. This proved to be sufficient to separate interstellar 
and telluric CII4 lines and guarantees their reality. The 
spectra presented in this paper are an average of all the 
data. The data were reduced in a standard way, using 
IDL routines (Boogert, Blake, & Tielens 2002). Atmo- 
spheric absorption features were divided out using the 
standard star HR 8585 (AlV), which is bright (^=3.78) 
and reasonably close to NGC 7538 : IRS9 (an airmass 
of 1.20-1.30 versus 1.34 for NGC 7538 : IRS9). The 
spectral shape and hydrogen absorption features in the 
standard star were divided out with a Kurucz model at- 
mosphere. The hydrogen absorption profiles in the Ku- 
rucz model are accurate to 1-2%, and residuals can be 
seen near strong lines. This does not affect the region 
of main interest for this paper (3.15-3.65 /im), except 
perhaps near 3.297 /im due to the HI Pf S line. Over- 
all, a good telluric correction was achieved, resulting in 
signal-to-noise values of ~110 in the region of CH4 at 
3.32 /im, ~ 40 at the bottom of the H2O ice band at 
3.0 /im and ~ 200 at wavelengths > 3.5 /im. Strong 
atmospheric lines, i.e. with less than 50% of the max- 
imum transmission in each setting, leave residuals and 
were removed from the final spectrum. The spectra were 
wavelength calibrated on the atmospheric emission lines, 
using an atmospheric model spectrum (P. van der Valk 
and P. Roelfsema, priv. comm.). All orders were sub- 
sequently combined by applying relative multiplication 
factors. For the analysis of the emission lines fi )3.2.1|l . 
the spectrum is flux calibrated by scaling to the contin- 
uum brightness of the data presented in Whittct ct alj 
l(1996f) . Finally, in conjunction with the CH4 features, we 
analyze Keck/NIRSPEC M ba nd CO lines in NGC 7 538 
: IRS 9, previously published in lBoogert et alJ l)2002|l . 



2.2. ISO/SWS R=l, 500 Spectroscopy 

The 3.26-3.38 /im spectral region of NGC 7538 : 
IRS 9 was observed multiple times in the high resolu- 
tion R= 1,500 AOT06 mode of the Short Wavelength 
Spectromet er on board the Inf r ared Space Obser vatory 
(ISO/SWS; IKessler et al1[T9M |deGraauweraOE996). 
From the ISO archive we have selected the longest inte- 
grations in this wavelength range, corresponding to TDT 
numbers 39002336 (UT 11 Dec 1996), 75101049 (UT 05 
Dec 1997), and 85200455 (UT 16 Mar 1998). AU obser- 
vations were part of the ICE_B AND/DWHITTET guar- 
anteed time program, and have not been previously pub- 
lished. The combined on-target integration time in the 
3.26-3.38 /im region amounts to 5219 seconds. The in- 
dividual data sets were processed with SWS Interactive 
Analysis and calibration flies version 10.1 in April 2003 at 
SRON/Groningen. The three AAR products were com- 
bined and the detector scans were aligned to the mean 
using second order polynomials. Then, in this highly 
oversampled spectrum, data points deviating more than 
3(7 from the mean per resolution element were removed. 
Finally, the data points were averaged and re-sampled to 
an R=l,500 Nyquist sampled grid. The resulting S/N is 
100, but signiflcantly poorer at the outer 0.015 /im edges. 
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Fig. 1.— The R=25,000 Keck/NIRSPEC L band spectrum of NGC 7538 : IRS9 on an optical depth scale, after correction for H2O 
ice absorption. The broad, deep absorption is the well known, yet unidentified, "3.47 /im" absorption band. The presence of absorption 
by the C-H stretching mode of solid CH4 is confirmed in the R=l,500 ISO/SWS spectrum. P, Q, and R-branch wavelengths of gas phase 
CH4 arc indicated at the top, where the numbers refer to the lower rotational J values. The gas phase lines are further highlighted in 
Fig. 1^ Ad ditional solid sta te absorption features are superposed on the 3.47 fim band. The most secure identification is that of methanol 
^llamandola et al, 11199211). Other possible identifications are indicated by arrows, but these are not further discussed in this paper: 
ethanol lAUamandola ct al. 1992), ethane iBoudin et al. 1998), formaldehyde iBchutte et al. 1996), and diffuse medium aliphatic C-H 
stretching modes ^Sandford et aljil991iV The weak, broad inflection near 3.25 /im is possibly due to absorption by aromatic hydrocarbons 
IBrooke et alj|1999l) . 



3. RESULTS 



3.1. Features in the 3.1-3.6 /ini Spectral Region 

The L band spectrum of NGC 7538 : IRS9 is domi- 
nated by the deep a bsorption band of the 0-H stretch- 
ing m ode of H2O ice llWillner et al.lll982HWhittet et al. I 
119961) . Its peak optical depth at 3.0 /xm is ~ 4.5. 
The most prominent substructure in the long wave- 
length wing of H2O is the well known, yet unidenti- 
fied, 3.47 /xm absorption band iAllamandola et al. Il992t 
IBrooke. Selkren. fc SmithI Il996t iDartois et all 1200^ 
We applied a smooth spline continuum to extract this 
band on an optical depth scale (Fig.^l. The results were 
verified by applying low order polynomials to contin- 
uum regions described in previous studies (e.g. Brooke, 
Sellgren, & Geballe 1999). The absorption shape is 
in good agree ment with th e low resolution spectrum of 
lAUamandola c t al. I l|1992() . as is the presence of a sub- 
band at 3.53 /im, attributed to the C-H stretching mode 
of solid CH3OH. A previously undetected absorption 
band is apparent by a downturn in the continuum slope 



at 3.32 /im. Although it is in a spectral region badly 
affected by the telluric CH4 Q-branch, its presence is 
confirmed in Keck/NIRSPEC observations on different 
dates. In addition, a feature at the same wavelength is 
found in the space based ISO/SWS spectrum (Fig. 
We thus conclude that this feature is real and is due to 
the stretching mode of solid CII4 fi l3.3|l . Not immedi- 
ately obvious in Fig. ^ but clearly detected upon close 
inspection (Fig. [2J) are P, Q, and R-branch lines of gas 
phase CH4. It is the first time that the 3.32 /im ice and 
gas features are detected in a circum-protostellar envi- 
ronment. 

A number of additional features, not further dis- 
cussed in this paper, are visible in the high resolution 
Keck/NIRSPEC spectrum. Although they are repro- 
duced in observations at different dates, they are en- 
tirely new and need to be confirmed in a larger sample 
of sources, or, as for solid CII4, in independent data. A 
weak feature at 3. 36 /im coincides with the C-H stretch 
modes of ethanol ijAllamandola et al. 11199^ and ethane 
(Boudin, Schutte, & Greenberg 1998). Interestingly, an- 
other band of ethane at 3.47 /im seen in laboratory ices 
appears to be present in NGC 7538 : IRS9 as well. Fi- 
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Fig. 2. — Magnification of the Keck/NIRSPEC spectrum of NGC 7538 : IRS9 in spectral regions surrounding gas phase CH4 ro- vibrational 
lines (AA=0.006 /im or /\v= 530 km s-i). This is the same data as shown in Fig. however on a flux/continuum scale. The gaps in 
the spectra are regions with less than 50% atmospheric transmission that were removed from the data. The thick gray line superposed 
on each observed spectrum depicts a two component LTE model, composed of a 6 km s~^ blue-shifted absorption component (T = 55 K, 
N = 2.5 X 10^® cm~2), and an 8 km red-shifted emission component (T = 600 K, N = 5.8 X 10^^ cm"^; Table|l}. The ro-vibrational 
transition, the lower rotational energy level (in K), and the center wavelength (in ^tm) are indicated in the lower left and right corners of 
each panel respectively. Note that CH4 emission is detected up to the highest excitation levels, i.e. up to R(15) and P(14). 



nally, weak absorptions at 3.47 and 3.41 /xm are also 
close to the C-H stretch mode of single bonded carbon in 
CH and CH2 groups in aliphatic hydroc arbons, normally 
seen in diffuse lines of sight only (e.g. iPendleton et al.l 
ir994) . 

3.2. Gas Phase Line Profiles at R^25,000 



3.2.1. Gas Phase GH4 

Ro-vibrational absorption and emission lines of gaseous 
CH4 are clearly detected toward NGC 7538 : IRS9 
(Fig. O. Although the telluric CH4 lines are deep, the 
large blue-shift of NGC 7538 : IRS9 with respect to the 
earth (§ 12.1(1 separates many interstellar lines well. The 
lines appear mostly in absorption up to lower energy lev- 
els with J = 3, and entirely in emission above J — 6. 
Intermediate J levels have P Cygni-likc line profiles, i.e. 
red-shifted emission and blue-shifted absorption. Such 
line profiles are usually associated with expanding en- 
velopes (see i)4.2l however). Note that for J > 4 multiple 



emission and absorption peaks appear for each J level. 
These line splittings are intrinsic to the CH4 molecule 
and result from rotation-vibration interactions with the 
spin manifolds of the four equivalent H nuclei. Thus E~ 
, F-, and A-typc multiplcts occur with spin statistical 
weights of gs=2, 3, and 5 respectively. 

The observed lines are analyzed using the 2001 u pdate 
of the HITRAN 2 000 database fRoth man eraP 120031: 
iBrown et aDl2003(l . All HITRAN "transitions contribut- 
ing to each observed interstellar line arc identified and 
the number of transitions responsible for at least 90% 
of the line strength is indicated in column 3 of Tables 
1 and 12 For example, at the resolution of our observa- 
tions, the P(ll) transition is split in 4 separate emission 
lines, which are composed of some 9 strong transitions. 
The contribution of each observed line to the total emis- 
sion or absorption per J level is indicated in column 4 
of Tables 1 and [51 and is calculated in the optically thin 
limit as 



9nJ-^nJ 



Detection of the Methane Stretch Mode 



TABLE 1 

CH4 Emission Line Strengths toward NGC 7538 



IRS 9 



Av{J) 




lines'' 


% 


10-21 w/cm 


R(15) 


3.15812-3.15864 


11 


84 


5.2 (1.1) 


3.15767-3.15792 


3 


13 


<3.3 


R(14) 


3.16677-3.16720 


13 


73 


3.8 (0.7) 




3.16626-3.16665 


5 


26 


1.5 (0.4) 


R(13) 


3.17545-3.17595 


13 


81 


9.0 (0.8) 


3.17498-3.17534 


4 


16 


4.3 (0.8) 


R(12) 


3.18435-3.18485 


14 


74 


13.6 (3.0) 




3.18390-3.18428 


4 


26 


4.8 (0.5) 


R(ll) 


3.19295-3.19325 


2 


21 


3.8 (1.5) 


R(10) 


3.20211-3.20253 


3 


34 


3.5 (1.0) 


R(9) 


3.21178-3.21229 


6 


78 


10.7 (3.0) 




3.21146-3.21179 


2 


22 


3.0 (0.8) 


R(8) 


3.22087-3.22118 


3 


47 


4.8 (0.8) 


R(6) 


3.24014-3.24080 


6 


100 


25.1 (1.5) 


R(5) 


3.25019-3.25064 


4 


100 


12.5 (1.4) 


R(4) 


3.26017-3.26054 


4 


100 


7.5 (1.0) 


R(3) 


3.27035-3.27063 


3 


100 


6.6 (0.5) 


R(l) 


3.29105-3.29127 


1 


100 


1.5 (0.5) 


R(0) 


3.30165-3.30187 


1 


100 


1.5 (0.6) 


P(2) 


3.33435-3.33463 


2 


100 


1.8 (0.4) 


P(3) 


3.34573-3.34600 


3 


45 


2.1 (0.5) 


P(4) 


3.35708-3.35736 


1 


23 


1.9 (0.6) 


3.35679-3.35703 


3 


77 


1.9 (0.9) 


P(5) 


3.36860-3.36901 


2 


45 


>4.8 


3.36814-3.36860 


2 


55 


4.5 (0.6) 


P(6) 


3.37967-3.38022 


3 


48 


7.8 (0.8) 


P(7) 


3.39182-3.39230 


2 


42 


> 5 




3.39138-3.39176 


2 


32 


3.5 (0.4) 


P(8) 


3.40379-3.40411 


1 


14 


>1.6 


3.40318-3.40370 


3 


48 


6.8 (0.4) 


P(9) 


3.41577-3.41613 


2 


18 


> 2.7 


3.41509-3.41560 


2 


22 


4.1 (0.4) 


P(10) 


3.42719-3.42772 


3 


35 


4.2 (0.7) 


P(ll) 


3.44128-3.44172 


3 


30 


3.2 (0.4) 


3.44072-3.44109 


1 


11 


1.6 (0.2) 




3.44009-3.44057 


3 


37 


3.9 (0.2) 




3.43944-3.43982 


2 


22 


2.0 (0.2) 


P(12) 


3.45380-3.45441 


4 


43 


4.7 (0.3) 




3.45313-3.45347 


2 


13 


1.4 (0.3) 




3.45252-3.45296 


2 


16 


1.1 (0.4) 




3.45175-3.45208 


3 


27 


1.9 (0.3) 


P(13) 


3.46664-3.46711 


3 


23 


1.1 (0.4) 




3.46566-3.46629 


3 


31 


0.8 (0.3) 




3.46501-3.46548 


3 


28 


1.8 (0.3) 




3.46406-3.46457 


2 


17 


0.8 (0.3) 


P(14) 


3.47952-3.48004 


3 


22 


1.9 (0.4) 




3.47814-3.47899 


4 


35 


0.6 (0.3) 




3.47757-3.47796 


2 


16 


< 0.8 




3.47659-3.47707 


3 


27 


1.7 (0.4) 



^ integration region after correction for 8 km s ^ red- 
shift 

number of blended lines in integration interval 
percentage of emission with respect to total emission 
for this upper J level 
Icr errors in brackets 



where gij is the product of the rotation (2J+1) and nu- 
clear spin (gs) statistical weights and Aij is the corre- 
sponding Einstein coefhcient for transition i at rotational 
level J . The summation over i is for all transitions in the 
wavelength range of an observed emission line (column 
2 in Tables 1 and[2Jl while the summation over n is for 
all transitions at rotational level J. For several J levels 
only part of the emission or absorption is observed as 
some components lie in poor atmospheric regions. For 
example, only half the R(8) line flux is seen and for 
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Fig. 3. — Rotation diagram of CH4 emission lines toward 
NGC 7538 : IRS9. P-branch lines are represented by filled cir- 
cles and R-branch lines by open circles. The solid line repre- 
sents a least-squares best fit to the data corresponding to a ro- 
tational temperature of 600 K at a column density of 5.8 xlO''^ 
cm~2 in a solid angle of 4.3x10— sterad. The dashed lines 
represent extreme (3(t) physical parameters that still fit the data: 
Trot =500 K and 7V(CH4)=20x 10^^ cm'^, and Trot =700 K and 
Ar(CH4)=2.6xl0i^ cm-2. 



the Q-branch the clearly detected absorption component 
represents only a few percent of the total Q-branch ab- 
sorption. Finally, the observed emission line fluxes Fj 
listed in Table 1 are corrected for the underlying H2O 
absorption using a w ide coverage ISO/SWS spectrum 
ijWhittet et al. iri996(l . The resulting continuum flux of 
1.5 Jy is corrected for foreground extinction using the 
Galactic extinction law IjDra inc 2005) and the Ay de- 
rived in ti4.2l Thus, the observed line fluxes are scaled 
to an uncxtincted continuum flux of 23.7 Jy. 

In order to derive basic physical parameters, we con- 
struct a rotation diagram for the CH4 emission lines. 
The Boltzmann equation of a single transition i at upper 
rotational level J is 



9iJ 



(2) 



where Nij and Eij are the column density and the energy 
of the upper level of a single multiplet component respec- 
tively. Assuming that coUisional excitation is applicable, 
Eij inc ludes a vib rational v = 1 excitation energy of 
4348 K l)Liddll994j) . TVtot and T^ot are the total CH4 col- 
umn density and rotational temperature. Throughout 
this paper we take the ro-vibrational partition function 
Q{Trot) from the HITRA N database, which i s applicable 
at Trot = 70 - 3000 K l|Fischer et all 12003(1 . At lower 
temperatures we calculate Q(Trot) from lines in the HT 
TRAN database. In the optically thin limit, the flux F'j 
of a partial J transition is obtained by summation over 
all multiplet lines 



i 

and thus the observed Fj relates to eq. [3]as 



(3) 
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TABLE 2 

CH4 Absorption Line Strengths toward NGC 7538 : IRS9 



Av{J) 


A'' 


lines'' 


% 


t>^, WlV_i.Lrll 

10-3 cm-1 


R(4) 


3.26003-3.26031 


3 


77 


10.0 (1.4) 


R(3) 


3.27019-3.27051 


3 


100 


12 (1) 


R(l) 


3.29092-3.29121 


1 


100 


12.7 (0.4) 


R(0) 


3.30151-3.30179 


1 


100 


17.6 (1.0) 


Q 


3.31234-3.31296 


4 


2 


35 (6) 


P(2) 


3.33423-3.33453 


2 


100 


6.5 (0.8) 


P(3) 


3.34530-3.34584 


3 


100 


12.3 (0.7) 


P(4) 


3.35660-3.35694 


3 


61 


5.9 (0.8) 


^ integration region after correction for f 
shift 


) km s— ^ blue- 



^ number of blended lines in integration interval 

percentage of absorption with respect to total ab- 
sorption for this lower J level 
la errors in brackets 



^ Mot ^-Ej/kT,„, /^x 

hCKj xnxJ2,9^JA^J Q(T,ot) ' ^ ' 

Here we make use of the fact that the upper energy level 
Ej and transition frequency vj are nearly the same for all 
multiplet lines (to within 0.1%). 11 is the size of the emit- 
ting region, for which we assume the value of 4.3x10"^'' 
sterad derived from CO emission lines, two orders of mag- 
nitude smaller than the solid angle corresponding to the 
0.42" wide slit (§|l!21l- The rotation dia gram constructed 
from the observed CH4 emission lines indicates a rota- 
tional temperature of 600±100 K (Fig. EJ. Ahhough 
Trot is well constrained, the derived column density of 
5.8^3^ X lO^^cm"^ is accurate to within a factor of ^3 
only (3(7 errors). 

To further illustrate the goodness of fit, to investigate 
possible optical depth effects and to model the complex 
blend of CH4 absorption and emission lines, we have also 
constructed single temperature LTE models. The level 
populations are calculated according to eq. [21 Emission 
line fluxes and absorption equivalent widths at any op- 
tical depth ar e calculated following Gcnzcl (1992), and 
iSpitzed 1)1978(1 respectively, applying the Voigt function 
for the line profiles. The observed line profiles are fit- 
ted by adding a single temperature red-shifted emission 
model to a single temperature, blue-shifted, absorption 
model. We find that a blue-shift of 6 km s~^ and a red- 
shift of 8 km s~^ are needed, and that both components 
are unresolved at the spectral resolution of 12 km s~^. 
The validity of the temperature and column density de- 
rived from the emission line rotation diagram is con- 
firmed by the relatively low model line optical depth; 
for example, at a Doppler parameter of 6d = 1 km s~^, 
Ntot = 5.8 X lO^*^ cm-2, and T=600 K, the peak optical 
depth is 0.33 (in the R(8) line). The absorption col- 
umn is well constrained to iV(CH4)=3±l x 10^^ cm~^ 
at Tex=55±15 K for intrinsic line widths of 6d=0.5- 
3 km (Table EJ. The narrowest lines are optically 
thick (t=5 at 6d=0.5 km s^^), while they are opti- 
cally thin for 6d> 1 km s^^. We take 6d=3 km 



as a preferred value throughout this paper, because the 
larger line widths are us ually p referred toward NGC 7538 
: IRS 9 l|Mitchell et al. Ill990t fBoonman fc van Dishoeckl 
120031) . but values of &d> 3 km s~^ are inconsistent with 
the unresolved CH4 lines observed. The physical param- 
eters derived from the CH4 stretching mode are in good 
agreement with the parameters from t he bending mode 
at 7.67 jim in lower resolution spectra l)Lacv et alJll994 
IBoogert et all lTgM). 

Finally, absorption or emission lines of the ^^CH4 and 
CH3D isotopes have not been detected. Study of the 
^•^CH4 lines is complicated, because the isotope wave- 
length shift is comparable to the separation of J levels 
and thus the ^^CH4 and ^^CH4 lines are blended. The 
column density upper limits are a factor of 3-7 smaller 
compared to the CH4 column density. At ^'^C/'^^C 
and H/D isotope ratios of at least 50 in de nse clouds 
l|Wilson fc Rood 1119941 ICeccareUi et al.ll200H) . these up- 
per limits are not significant. 

3.2.2. Gas Phase CO 

At first glance, the line profiles of the fundamental 
ro-vi brational CO trans ition s at 4.7 /xm, as published 
in Mi tchell et al. 1 1| 199(1) and IBoogert et a l. (2002), look 
similar to those of CH4. They show P Cygni profiles as 



TABLE 3 

Gas Phase CO Absorption toward NGC 7538 : IRS9 





eq. 


width'" 




10" 


3 cm"! 




13 CO 


CI8O 


R(18) 


2 (1) 




R(17) 


5 (2) 




R(16) 


3 (1.1) 




R(15) 


6.5 (2) 




R(13) 


10 (2) 




R(12) 


16 (3) 




R(ll) 


12 (4) 


<4.5 


R(10) 


16 (3) 


<4.5 


R(9) 


18 (5) 


5 (1) 


R(7) 


30 (8) 


2.2 (1.0) 


R(6) 


30 (2) 


3.9 (1.2) 


R(5) 


41 (3) 


4.7 (1.2) 


R(4) 


71 (4) 


5.6 (1.0) 


R(3) 


71 (4) 


5.2 (0.8) 


R(2) 




13.7 (1) 


R(0) 


55 (3) 


10 (2) 


P(l) 


55 (2) 


6.4 (1.5) 


P(2) 


56 (3) 




P(3) 


74 (2) 




P(4) 


27 (2) 




P(18) 


<15 




P(19) 


3.4 (0.4) 




P(21) 


2.5 (0.8) 




P(22) 


< 6 





^ transitions not listed have 
not been observed, are in a 
poor atmospheric region, or are 
blended with a brighter isotope 
^ la errors in brackets 
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Fig. 4. — Line profiles of CH4 and CO isotopologues on a velocity scale corrected for the earth and heliocentric velocity. Lines of similar 
lower energy level are plotted: K in panel (a), and ~100 K in panel (b). For reference we have indicated the systemic velocity with a 
vertical dashed line. Note the very different profile of CH4 compared to CO lines. The -16 km s~^ outfiow is absent in CH4, an d th e 
R(0) CH4 line does not peak at the systemic velocity, indicating that CO and CH4 have different a chemical origin and evolution ( i|4.1i . 
Panel (c) shows emission-dominated CH4 and CO line profiles. The high- J CO line profile toward NGC 7538 : IRS9 is an average over the 
individual lin es P (27), P(28), P(30), and P(34). The emission lines clearly peak at -|-8 km opposite to what is expected for P Cygni 
line profiles ( M.2\ . 




200 400 600 800 1000 1200 
Ej/k (K) 

Fig. 5. — Rotation diagram of ^^CO (filled circles) and C^*0 
(open circles) absorption fines toward NGC 7538 : IRS9. The C^^^O 
column densities have been scaled u p by a factor of 6.9 to match 
the abundance of interstellar ^^CO 1 Wilson fc Ro od IIH99I). The 
good agreement of these isotopes in this plot implies that the ^^CO 
absorption lines are optically thin. The solid line represents a two 
component gas with rotational temperatures of 20 and 270 K both 
at a I^CO column density of 4.65 xlO^'' cm-^. Note that the 20 
K gas absorbs at the systemic velocity, while the warm component 
absorbs at much higher velocity ("Fig. El. 



well, in ^^CO, ^^CO, and perhaps even the C^^O iso- 
tope, with the strength of emission versus absorption 
increasing for higher J levels. Upon closer inspection, 
however, the absorption profiles of CO and CH4 differ. 
The low- J isotopic CO lines peak at km s~^, compared 
to -6 km s-i for CH4 (Fig. Hi). Higher- J CO hnes 
(ii^iow > 100 K), on the other hand, are dominated by 
absorption at blue-shifted velocities > 16 km s~^, which 
is absent in CH4 (Fig.^jD). This provides clues to the for- 
mation and evolution history of CH4 fi) l4.1|l . and in this 
context it is relevant to quantify the physical conditions 
of the absorbing and emitting CO gas as well (Table 0J. 

We analyzed th e R=25,000 Keck/NIRSPEC spectrum 
of iBoogert et al.l (^002) and arrive at different conclu- 
sions compared to Mitchell et al. (1990, 1991). First, 
and most importantly, in the Keck/NIRSPEC data blue- 
shifted absorption at 16 km s~^ is seen up to J-levels 
as hig h as 20, but only up to J = 3 in Mitchell et a^ 
(jl99ll). In the rotation diagram constructed from the 
integrated ^^CO absorption equivalent widths, this high 
velocity gas is responsible for the linear part above ener- 
gies of 100 K (Fig.[5j Table|21). From this a column den- 
sity of iV(i2co)=3.5±1.0 X 10^^ cm-2 at a temperature 
of T„t = 270 ± 100 K (3cr) is deri ved. Th e column den- 
sity is comparable to that found bv lMitchdl et al. ( 1991), 
but Trot is a factor of 10 higher in the Keck/NIRSPEC 
data. A second difference is that C-'^^O lines are detected 
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in the high S/N Keck/NIRSPEC data. The i^cO/C^^O 
absorption Hne ratios of corresponding J levels are equal 
to the interstellar isot ope ratio of 6. 9 at large Galacto- 
centric radii fWilson fc Rood I Il994|) within a 20% er- 
ror margin (see the C^^O data points in the ^^CO ro- 
tation diagram; Fig. 1^1 ■ An important consequence is 
that reliable CO column densities can be derived, and 
thus reliable CH4 abundances in the various velocity 
ranges (§4.1). The steep part of the ^^CO rotation di- 
agram, corresponding to gas at velocities of -8 to -f 11 
km s~^, indicates T^t =20±7 K and a column density 
of iV(i2cO)=3.5±1.0 X IOI8 cm -^ which is a factor o f 
4 less compared to that found bv lMitchell et al. 1 1)199(11) . 
This is due to our factor of 2 smaller measured equivalent 
widths and, presumably, lower assumed line opacity. 

For the emission lines, the observed ^^CO/^'^CO peak 
strength ratio is 6 ± 2 at intermediate upper J levels of 
7-15. Using the standard radiative transfer equation, the 
peak optical depth of the ^^CO emission lines, T12 can be 
derived from 

f (^^CO) _ 1 - e-^i^ 

F(i3C0) ^ 1 -e-^12/80' ^ ' 

At a i^co/^^co isotope ratio of 80 l|Boogert et all 
l2002j) . it follows that T12 = 13 ±4. Such an optical depth 
corresponds to a column density of A^(^^CO)=3.2± 1.0 x 
10^^ cm~^ at an intrinsic line width of 6d=3 km s~^. We 
further discuss the origin of the emitting gas in ii4.2l 

3.3. The Solid CH4 Absorption Bands 

Solid CH4 has been detected by its C-H bending mode 
at 7.67 um in ISO/SW S observations of NGC 7538 : IRS9 
(jBooeert et al.llT996() . The peak position and width of 
the interstellar band are well fitted with an amorphous, 
low temperature (10 K) laboratory ice in which CH4 is di- 
luted in an H2O matrix. This same laboratory ice gives 
an excellent fit to the new 3.32 /xm absorption feature 
detected in the Keck/NIRSPEC and ISO/SWS spectra, 
not only in peak position and width, but also to the 
peak optical depth of 0.1. We thus attribute this inter- 
stellar band in its entirety to the C-H stretching mode of 
solid CH4 in an Il2 0~rich ice (Fig.EJ. A mismatch with 
the observed ISO/SWS spectrum on the low wavelength 
wing is attributed to emission by the strong Q-branch 
of gaseous CH4 (Fig. 0) . Indeed, a better match is ob- 
tained after subtracting the model constructed from the 
observed P and R-branch hues from the ISO/SWS spec- 
trum (§3.2.1). In contrast, contamination of gas phase 
emission to the 7.67 /im ice band is small because the 
continuum emission at this wavelength is more than an 
order of ma gnitude brighter compared to the continuum 
at 3.32 /im l|Whittet et al. l ll996'). Indeed, a weak inflec- 
tion at 7.66 fim is due to Q-branch absorption. Simulta- 
neous fits of the C-H stretch and bend modes highlight 
the good fit of CH4:H2 ices (Fig. Pure CH4 ices and 
mixtures with CO, otherwise highly abundant in th is line 
of sight l|Tielens et al. Ill991t iBoogert et alJl2002D . can 
be excluded at a high confidence level. Note that the C- 
H stretching mode of CH4 better discriminates between 
polar and apolar ices than docs the bending mode. Mix- 
tures with other high-dipole molecules, such as NH3 and 
CH3OH, do not provide equally good fits compared to 




3.28 3.30 3.32 3.34 3.36 7.50 7.55 7.80 7.65 7.70 7.75 7.80 



Fig. 6. — The stretching mode of sohd CH4 observed with 
ISO/SWS (thick histogram) and KECK/NIRSPEC (thin his- 
togram), compared to a number of laboratory ice simulations in 
panel (a). At the very top the effect of removing the Q-branch 
emission of gas phase CH4 is shown. The CH4 bending mode in 
panel (b) is fitted with the same laboratory ice, without scaling 
the peak optical depth. Best fits are obtained with H2 0-rich ices. 
Note the poor fit of CO-rich ices to the CH4 stretching mode. 



H2O, but at present we cannot exclude combinations of 
these polar species. 

4. DISCUSSION 



4.1. Astrochemical Implications 

The detections of both the 3.32 /im stretch and the 
7.67 /im bend modes firmly establish the presence of 
both gaseous and solid CH4 in the interstellar medium. 
The gas-|-solid state abundance averaged along the line 
of sight of the massive protostar NGC 7538 : IRS9 is 
X[CH4]=2 X 10~^, and the absorption gas/solid state 
abundance ratio is 0.23. A combined high resolution 
CO/CH4 line profile analysis shows that the abundance 
as well as the gas/solid state abundance ratio of CH4 
vary locally by an order of magnitude, however (Table^. 
An ice-|-gas abundance of a few 10"^ is found in the 
outer envelope (systemic velocity), the inner envelope 
(-6 km s~^) and the compact emission component (-1-8 
km s^^). In the outer envelope, at 20 K the coldest com- 
ponent along the line of sight, at least 90% of the CH4 
is frozen on grains. Within the envelope a segregation is 
present as well; CH4 is intimately mixed with H2O and 
CH3OH rather than with CO. In the warm, low velocity 
inner regions of NGC 7538 : IRS9, CH4 is mostly present 
in the gas phase. These enhanced abundances strongly 
contrast to the absence of CH4 in the massive, high ve- 
locity -16 km s^^, warm (T ~ 270 K) outflow seen in 
CO (X[CH4]< 0.4 X 10^^). 

These new results provide clues to the formation and 
evolutionary history of CH4, refining previous scenarios 
(Lacy ct al. 1991; Boogert et al. 1998). In time depen- 
dent pure gas phas e models, such as the 'standard' dense 
cloud model of .Lee. Bettens. fc Herbsti fl996i) a nd th e 
UMIST model of iMillar. Farauhar. fc Willacvl l(T997|l . 



Detection of the Methane Stretch Mode 



9 



TABLE 4 

CH4 AND CO Physical Parameters toward NGC 7538 : IRS9 



Species Location'' 


r^HcLSe 


V 


cx 


iV 


/\r/ /\/"^^d,c 
JV/iVH 


IN OLGS 








km s ^ 


It^ 
Iv 


lu cm 


iu 




CH4 


outflow 


gas/abs. 


— 16 


270 


< 0.1 


< 0.4 


assumed T 


CH4 


core 


gas/abs. 


-6 


55±15 


0.3±0.1 


> 3 


f)D=0.5-3 km s-i 


CH4 


envelope 


gas/abs. 





20 


< 0.1 


< 0.3 


assumed T 


CH4 


core 


gas/em. 


+8 


600±100 






4.3 X 10-1'' sterf 


CH| 


envelope 


ice/abs. 






1.3±0.1 


2.1±0.7'' 




CO 


outflow 


gas/abs. 


-16 


270±100 


35±10 


100 


opt. thin i^CO 


CO 


core 


gas/abs. 


-6 


55 


< 10 


100 


assumed T 


CO 


envelope 


gas/abs. 





20±7 


35±10 


100 


opt. thin "CO 


CO 


core 


gas /em. 


+8 


500±150 


30±10' 


100 


4.3 X 10-14 sterf 


COJ 


envelope 


ice/abs. 






18±2 


36±10'' 





^ assumed location from velocitv. see also t) il4.ll and l4.2l ^ in gas or ices, in emission or absorp- 
tion; velocity w.r.t. systemic Vis,-=57 km s~^; errors are 3a; ° assuming A'^h/CO=1x10* 



from IT^cv et all 119941 ' solid angle derived from CO emission; B IBoQgert eTalHTogg: 
ice abundance with respect to CO gas+ice abundance in envel ope: ' uncertainty in column 
dominated by uncertainty in line opacity; J IBoogert et al.|[2002l 



gas phase CH4 abundances of a few 10"^ are achieved 
after 10^ years. Although this is comparable to the ob- 
served abundances, there are strong arguments against 
a gas phase origin: the low gas/solid state ratio and 
the intimate mixture between CH4 and H2O in the 
ices in the envelope. For comparison, CO, a typical 
gas phase product, has an order of magnitude larger 
gas/solid phase ratio compared to CH4, even in the cold 
gas phase. H2O on the other hand is formed by hydro- 
gena tion of oxygen on the g rains at low temperatures 
(e.g. I.Tones fc W illiam^'1984^. and a low gas /solid state 
ratio is observed (Boonman & van Dishoeck 2003). It is 
thus likely that CH4 was formed on grain surfaces by hy- 
drogenation of atomic C, at the same time as H2O was 
formed from atomic O. The low observed CH4 abundance 
points to a formation time whe n most of th e C was al- 
ready locked up in CO (Hoog-er t'^nnrmfll . 

The gas phase CH4 abundance in the -6 and -1-8 
km s-^ components is comparable to the solid state 
abundance in the envelope. The CH4 may thus have 
been released from the grains in the inner regions near 
the protostar. Because CH4 is diluted in an H20-rich 
ice (Fig.O, its sublimation temperature is close t o that 
of H2 O, 90 K [Tielens fc Whittet 1996,; Hirao ka et all 
^23)- The temperature of the +8 km emission com- 
ponent is indeed well above 90 K. The gas at -6 km s"^, 
however, has a temperature of 55±15 K only. While in 
one massive p rotostar, W 33A, C H4 gas at T = 110 K 
was detected l|Boogert et al.llT99^ . a lower temperature 
of 50 K was also measured in another line of sight, GL 
7009S (|Dartois et al.lll99a() . We speculate that the icy 
mantles did not thermally evaporate as a result of pro- 
tostellar radiative heating, but by a mild shock that has 
briefly heated the grains above 90 K. Subsequently the 
gas has cooled radiatively or adiabatically as the inner 
NGC 7538 : IRS9 core expands outwards with a (cur- 
rent) velocity of 6 km s"^. 

Shocks may have had a profound impact on the higher 
velocity gas. CH4 is absent at velocities of 16 km s^^ and 
above, where warm CO (T = 280 K) is very abundant. 
After release from the grains, CH4 must have efficiently 



chemically reacted. The observed CO velocities up to 
35 km s^^ point to the presence of non-destructive C- 
type shocks. These shocks have been well studied in the 
context of molecular chemistry. At initial densities of 
10^ cm^"^, shocks with velocities of 15-35 km s""'^ have 
peak temperatures of 500-2500 K (Bergin, Neufeld, & 
Melnick 1998). At such high temperatures, and at the 
low expected sh ock pressures'^, CH4 is rapidly driven into 
CO (jPrinn fc Feglevlll989^1 . The occurrence of a C shock 
therefore explains the absence of CH4 in the outflow 
of NGC 7538 : IRS9. The importance of shock chem- 
istry in high mass cores was dedu ced from low gas/solid 
CO2 abundance ratios as well ijCharnlev fc Kaufmanl 
l200^^■ and by the abu ndance ratios of S bearing species 
ijHatchell fc Viti| l2002fl . One distinct characteristic of 
shock chemistry is the forma tion of H2 O at temperatures 
above 250 K (e.g. Bergin ( ^1111119981) . H2O columns 
comparable to CO are expected, i.e. 2.8x10^* cm^^ for 
the NGC 7538 : IRS9 outflow (Table l^. This is a fac- 
tor of 3.5 larger than the upper limit to the gas phase 
H^O abundance derived from ISO/SWS 6 fim spectra 
(^oonman & van Dishoeck 2003). High resolution 6 ^ui 
spectra are required to further investigate the presence 
of H2O in the line of sight of NGC 7538 : IRS9, the effect 
of blended emission and absorption components and the 
relation to shocks and outflows. 

4.2. Nature of the Low Velocity Warm Gas 

The elevated temperature of the low velocity, 6-8 
km s"^, absorption and emission components seen both 
in CO and CH4 locates this gas in the inner regions of 
the protostellar core. The present observations further 
constrain the origin of this material. 

The area of the radiating physical body can be cal- 
culated by matching the ^^CO column density derived 
from the isotope ratio (§ 13.2.2(1 to the emission line 
fluxes. The measured "'^^CO line fluxes need to be cor- 

^ logio{P) = — 13 bar assuming a factor of 2 density enhance- 
ment at a maximum temperature of 2500 K fo r a 35 km s-^ shock 
IKaufman fc NeufeldllT99l IBergin et al.|[T99l) 
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reeled for continuum absorption by foreground mate- 
rial. We derive a 4.7 /im continuum extinction mag- 
nitude of A4.7 = 1.8 ± 0.6 using the interstellar ex- 
tinction law of iDrain^ l)2003(l and an optical extinction 
Ay of 60±20. The la tter is derive d from the Ay/Nii 
relation of .Bohhn. Savage, fc Drak3 l|1978.) and the av- 
erage line of sight hydrogen column density from this 
work (Table El a nd that from the silicate absorption 
band ( Ticlcns ct al . Il99lh . The extinction corrected line 
fluxes are thus matched at a solid angle of 4.3x10^^'* 
stera d. At the distance to NGC 7538 : IR S9 of 2.8 
kpc ijCramnton. Georgelin. fc Georgelin II1978D this cor- 
responds to a radius of 70 AU. The similarity of CO 
and CH4 velocities and emission line widths (Fig. ^) 
indicates that both molecules are present in this same 
volume. 

The CO and CH4 gas detected on these scales may 
well be related to the dense material that was inferred 
on scales of < 300 — 600 AU from u nresolved 3 millime- 
ter continuum interferometer images (jvan der Tak et al. I 
iOOO). This intense emission is unrelated to the envelope, 
and instead may originate from a dense shell or circum- 
stellar disk. The P Cygni-like line profiles in both CO 
and CH4 may suggest that an expanding envelope is the 
preferred geometry. In these line profiles the absorption 
component is red-shifted by the shell expansion velocity, 
and the emission component is centered at the systemic 
velo city and is broadened by twice the expansion veloc- 
ity. iMitchell et al.l l)199lD argue that the CO emission 
component is blue-shifted due to the effect of deep ab- 
sorption by foreground material at the systemic veloc- 
ity. This hypothesis can be tested by the profiles of the 
highest J transitions, which should not have foreground 
contributions. We find that even the lines from the high- 
est rotational levels of CO are blue shifted by -1-8 km s~^ 
(Fig. Et). Furthermore, the CH4 emission lines peak at 
this same velocity and are also unresolved. Perhaps the 
NGC 7538 : IRS9 core moves within its envelope by 8 
km s^^ and the observed peak of the emission line is the 
center of the P Cygni profile. In this scenario the shell ex- 
pansion velocity is 14 km s~^ and an emission line width 
of 28 km s^^ is expected. Such broad lines are inconsis- 
tent with the observations as well. We conclude that the 
line profiles are not true P Cygni type profiles and thus 
are inconsistent with a spherically symmetric expanding 
envelope. It is unclear what the relation is between the 
55 K absorbing gas seen at —6 km and the 500 K 
emitting gas at -1-8 km s~^. More detailed dynamical 
models of circumstellar disks, accretion flows and out- 
flows are required. The possibility of photon exci tation 
needs to be considered as well. As noted in Mitchell et al.l 
(|199l5 the similarity of CO absorption and emission line 
profiles between a number of massive protostars suggests 
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that these lines trace a fundamental property of their in- 
ner cores. This is strengthened by the similarity of the 
CH4 line profiles. 

5. CONCLUSIONS AND FUTURE WORK 

Using high resolution L and M band spectroscopy we 
have shown that the gas phase CH4 abundance varies 
considerably in the line of sight of the massive protostar 
NGC 7538 : IRS9. This is a result of specific molecule 
formation and destruction processes: CH4 is formed on 
grain surfaces and subsequently evaporated from the 
grains in the inner regions of the protostar. Both evapo- 
ration by photon heating and by shocks likely have played 
roles. CH4 efficiently 'burns' to CO in the highest veloc- 
ity shocked regions. The gas phase CH4/CO abundance 
ratio is thus a sensitive tracer of shock chemistry. High 
spectral resolution observations of other species, in par- 
ticular H2O at 6 /im, are required to further test the 
importance and nature of shock chemistry in high mass 
star forming regions. The Stratospheric Observatory For 
Infrared Astronomy (SOFIA) will off er this opportunity 
in the near future ijLacv et alJl2002j) . Observations of 
the 3.32 /im CH4 stretch mode in more lines of sight 
require high radial velocities (> 30 km s~^) or temper- 
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rated telluric lines. In this respect, the 7.67 /im bending 
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the CH4/CO ratio is a sensitive tracer of protoplanets. 
In the relatively low pressure solar nebula low CII4/CO 
abundances occur. The higher pressure in protoplane- 
tary sub-nebulae, however, favors higher CII4/CO ratios 
fPrinn & Fcglcv 1989). Variations of the CII4 abundance 
between solar system comets possibly trace such physi- 
cal condition grad ients in the early protoplanetary disk 
lIGibb et alJl200l . 
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